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Abstract. - With reference to the broad and narrow Fano-Feshbach resonances of 6 Li at 
about 822G and 543G, we show that for the broad resonance a molecular coupled-channel 
calculation can be mapped with excellent accuracy onto an effective single-channel problem 
with a contact interaction. This occurs for a wide enough range of the magnetic field, that 
the full BCS-BEC crossover can be realized with a typical trap. For the narrow resonance, 
the mapping onto a single-channel problem and the realization of the BCS-BEC crossover are 
restricted to too narrow a range of the magnetic field to obtain them in practice. A general 
criterion is also formulated for deciding whether the BCS-BEC crossover can be exhausted 
within the single-channel model for any specific Fano-Feshbach resonance. In this way, the 
BCS-BEC crossover for Fermi atoms with the broad resonance is placed on the same footing 
as the corresponding crossover for different physical systems. 



Fano-Feshbach (FF) resonances [1] are currently used to control the effective atom-atom 
interaction in trapped Fermi gases [2-5], for realizing the BCS-BEC crossover [6-8] from 
overlapping Cooper pairs to non-overlapping composite bosons at low enough temperature. 
Study of this crossover was originally motivated by the condensation of excitons in solids [9] , 
and more recently applied to nuclei [10] and high-temperature superconductors [11]. 

In this context, a large amount of work has been made by adopting a contact potential for 
the effective fcrmion-fermion attraction, regularized in terms of the scattering length ap [12- 
14]. The use of this potential considerably simplifies the many-body diagrammatic structure, 
both in the normal [13, 14] and broken-symmetry [13, 15] phases. By this approach, only 
fermionic degrees of freedom are retained in the many-body Hamiltonian, in the same spirit 
of the original BCS theory [16]. Since just the quantity oj? is varied in a controlled way 
by sweeping the magnetic field across a FF resonance, it would appear that trapped Fermi 
gases constitute an ideal testing ground for many-body theories based on the above two-body 
interaction. Previous theoretical work using the same interaction could thus be adapted to 
trapped Fermi gases with limited effort. 

The use of the same effective two-body interaction for such different systems like high- 
tempcrature superconductors and trapped Fermi gases may also lead to the emergence of 
universal features. This should be especially desirable, as the insights gleaned from the 
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BCS-BEC crossover with trapped Fermi gases could lead to a better understanding of high- 
temperature superconductors [17]. 

Most theoretical work on the BCS-BEC crossover with trapped Fermi gases, however, has 
been formulated with a fermion-boson model, aiming at incorporating the molecular states 
coupled in a FF resonance [18]. Inclusion of resonance processes in the BCS-BEC crossover was 
then named "resonance superfluidity". This was claimed to result in a different BCS-BEC 
crossover [19] from that discussed for high-temperature superconductors. A fermion-boson 
model was actually proposed some time ago for the study of high-temperature superconductors 
[20]. This model was, however, conceived to include phenomenologically the coupling of 
fermions to a boson mode, which represents fluctuations internal to the fermion system (and 
not a two- fermion state as in resonance superfluidity). 

Given the current experimental advances on the BCS-BEC crossover with trapped Fermi 
atoms, it seems timely and important to settle the issue of which (single- vs multi-channel) 
model is relevant for an accurate description of this crossover. By performing ab initio calcu- 
lations of molecular 6 Li in the presence of a magnetic field, comparison will be made between 
the two FF resonances occurring at (about) 822G and 543G, which we shall regard as repre- 
sentative of "broad" and "narrow" resonances, respectively. We will show that, for the broad 
resonance, the outcomes of the coupled-channel calculation can be mapped with excellent 
accuracy onto the results of the single-channel model with a contact potential over a wide 
magnetic field range. On the contrary, for the narrow resonance this mapping can be realized 
only for too narrow a region of the magnetic field to be experimentally accessible. 

These two resonances have been selected to pursue the realization of the BCS-BEC crossover 
with trapped Fermi atoms (cf. Refs. [2] and [3,4] for the narrow and broad resonance, in the 
order). We will conclude that the experimental realization of the BCS-BEC crossover can 
be made, in practice, only with the broad resonance and that the use of the single-channel 
model is correspondingly appropriate to describe the crossover [21]. The ultimate reason for 
the "broad" and "narrow" nature of the FF resonances here considered will also be accounted 
for. 

The key theoretical issue is whether a minimal description of many-body effects in these 
systems can be given in terms of a single-channel model (where the scattering length ap is 
the only relevant parameter) or of a multi- (possibly two-) channel model (where additional 
parameters, like the effective range of the potential, appear). In the theory of resonance 
superfluidity [18], the need for a multi-channel model to describe the BCS-BEC crossover with 
trapped Fermi atoms was assumed, arguing generically that the energy dependence of two- 
body scattering matters for this problem [22]. The point is, however, that the experimental 
conditions introduce the average interparticle distance as an intrinsic length scale, so that 
only after combining this length scale with the two-body properties of a given FF resonance 
one can decide whether the single- or multi-channel model is adequate. 

The molecular calculation for the biatomic molecule of 6 Li is set up by the Born-Oppcnhcimcr 
scheme. The auxiliary electronic Schrodinger equation at fixed nuclear separation R is solved 
via a standard Configuration Interaction method, with zero projection of the total electron 
angular momentum along the nuclear axis. The coefficient of the R~ 6 tail of the electronic 
energy curves is fine-tuned (to within a few percent), in order to fix the positions of the FF 
resonances at the experimental values (822±3)G [3,4] and (543. 23± .05)G [2] (these positions 
turn, in fact, out to be quite sensitive to the details of the curves for the electronic energy). 
[In practice, the R~ 6 tail of these energy curves is set to zero for R larger than i? = 4000 
a.u..] With this slight adjustment only, the scattering length is found to have a zero crossing 
at B — 533G, which is quite close to the experimental value (528 ± 4)G [23]. In addition, the 
calculated hyperfme splitting between the atomic F = 3/2 and F = 1/2 states turns out to be 



S. SlMONUCCI, P. PlERIand G.C. STRINATi: BROAD VS NARROW FANO-FESHBACH RESONANCES3 




5 10 15 20 25 250 500 750 10001250 

R (a.u.) B(G) 

Fig. 1 Fig. 2 

Fig. 1 - Molecular energy curves plus (diagonal) hyperfine and Zeeman contributions (taken at 820G) 
for the singlet (S) and triplet (T) states of lowest energy vs the nuclear separation R. The inset 
shows the details of the crossing at R = 18 a.u.. 

Fig. 2 - Scattering length dj? vs magnetic field B. The inset shows the details of the narrow resonance 
at about 543G. 



217 MHz, to be compared with the experimental value 228 MHz. The nuclear wave equation 
contains the hyperfine and Zeeman couplings, and is solved only for s-wave relative motion. 
The low-energy sector relevant to the two FF resonances is thus spanned by five channels, 
corresponding to the singlet with spin configurations |0, 0) c |0, 0) n and |0, 0) c |2, 0) n (which 
will be referred to as the "resonance" channels), and to the triplet 3 S+ with spin configura- 
tions |1, l) c |l, l)n, |1,0) C |1, 0) n , and |1, l) c |l, l)n (which will be referred to as the "scattering" 
channels), for the electronic (e) and nuclear (n) spin functions with F z = 0. 

Figures 1 and 2 show representative outcomes of our molecular calculation. In particular, 
in Fig.l the (diagonal) hyperfine and Zeeman contributions (taken at 820G) are added to 
the electronic energy curves of lowest energy for the singlet (S) and triplet (T) . Note in the 
inset the crossing of the two curves induced by the Zeeman coupling, which makes the lowest 
energy threshold of triplet character. In Fig. 2 the calculated scattering length is shown 
from B — to B — 1300G. The inset shows the details of ap for the narrow resonance at 
about 543G. 

Several quantities of interest can be extracted from this calculation. Below resonance, 
we compare: (i) The binding energy Eh (obtained from the multi-channel calculation) to 
the value eo = (Map)^ 1 for the contact potential, where M is the nuclear mass; (ii) The 
corresponding (root-mean-square) radius R of the full bound wave function to the value ap/^/2 
for the contact potential. We also calculate: (iii) The (root-mean-square) radius R2 associated 
with the component of the total wave function in the resonance channel; (iv) The (squared) 
projection |u>i| 2 of the total wave function in the scattering channel; (v) The corresponding 
projection {(f> C p\wx) onto the bound wave function <p cp — (^/2/ap) exp(—R/ap) for the contact 
potential. 

Table I reports all these quantities for the two resonances. Note how, for the broad res- 
onance, the effective single- channel model with a contact potential reproduces with excellent 
accuracy all results obtained by the multi-channel calculation. Note, in particular, how the 
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Table I - Comparison of the molecular calculation with the effective single- channel model, for the 
broad resonance at about 822G and for the narrow resonance at about 543G. 
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bound-state wave function of the multi-channel calculation projects almost completely in the 
scattering channel, which is the sole considered by the effective single-channel model. The 
spatial extension R 2 of the boson introduced in resonance-superfluidity theory [18] remains 
instead quite smaller than the extension R of the true molecular wave function (which corre- 
sponds to the internal wave function of the composite boson). A different situation occurs for 
the narrow resonance. In this case, the effective single- channel model is appropriate only too 
close to the resonance for being of practical relevance. 

To study the BCS-BEC crossover with FF resonances, the above analysis of the molecular 
calculation must be complemented by the value of the Fermi wave vector fop (determined by 
the total number of Fermi atoms in the trap and the trap frequencies) and by the minimum 
experimental accuracy for the magnetic field. 

Irrespective of the underlying theoretical model, the dimensionless parameter (kpap)^ 1 
should exhaust the BCS-BEC crossover within a range 1 about the unitarity limit at 
(kpap)^ 1 = 0. This is because the BCS (with ap < 0) and BEC (with ap > 0) regimes 
should be reached when kp\ap\ <C 1, while in the crossover region kp\ap\ diverges. To span 
the crossover, one thus needs to identify (at least) three representative values (say, -1.0, 0.0, 
and 1.0) of (kpap) -1 , by tuning the magnetic field that controls the FF resonance. With kp = 
2-^3 x 10 -4 a.u. for the experiments of Refs. [3] and [4], these values of (A^of)" 1 correspond 
approximatively to the magnetic field values (1300, 822, 730)G, which are separated by a step 
SB ~ 100G much larger than the minimum experimental accuracy. For these values of B, the 
single-channel model is totally appropriate to describe the two-body scattering, as seen from 
Table I. The situation is reversed for the narrow resonance. In this case, the above values of 
(kpap)- 1 are realized by setting the magnetic field at (543.2225, 543.2217, 543. 2209)G, with 
a step SB ~ 0.001G fifty times smaller than the minimum experimental accuracy [2]. 

To complete the mapping into the single-channel model, there remains to identify the 
effective spin states | |) e ff and | |) e ff- From the molecular calculation, the wave function for 
the broad resonance is found to be (essentially) a triplet both for electrons and nuclei. This 
implies the identification | f ) c ff I5, — \)c 1 1 , 0} n and | J.) e ff <-» 1 5, — j) c |1, l) n for the separate 
atoms, as well as the effective singlet configuration (| f)^ | J,)^ — | f)^ | l)cff)/V / 2 for the 
pair of atoms (with the labels 1 and 2 referring to the two atoms). 

Further information can be extracted from the molecular calculation when the total energy 
is above threshold. It has been asserted [18] that, in resonance conditions, the scattering length 
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Fig. 3 - Mechanism for the occurrence of the narrow and broad FF resonances (see text). 

ap no longer suffices to account for the scattering properties, since these depend strongly on 
energy. Quite generally, the scattering amplitude takes the form f(k) = (g(k) — ik)~ l , where 
the wave vector k is associated with the kinetic energy above threshold. For a contact potential 
g(k) = —dp 1 , while for a multi-channel problem g(k) in the scattering channel contains 
additional fc-dependent terms. In particular, at the lowest order in k, g(k) = — aZ 1 + r^k 2 /2 
where r identifies the "effective range" of the potential. Near resonance (when a F ~ 0), f(k) 
has a strong fc-dependence even in the absence of ro. The term r$k 2 /2 begins to matter only 
when r k ~ 1. For the BCS-BEC crossover, a natural cutoff for k is provided by kp in the 
weak-to-intermediate coupling regime and by aj, 1 in the strong-coupling regime. The term 
r k 2 /2 is thus irrelevant provided |r |fc_F <C 1 and \r \ <§C |ap|. The values of r are reported 
in Table I for the two resonances. For the broad resonance, ro is positive (as it is the case for 
a potential problem) and remains much smaller than |of|. For the narrow resonance, none 
of these properties are verified. Since for the broad resonance the product r$kp is smaller 
than 1CP 2 (with a typical value kp as 10~ 4 a.u.), we conclude that the energy dependence of 
the scattering properties (over and above that resulting from a contact potential) is actually 
irrelevant when realizing the BCS-BEC crossover with the broad resonance. 

It is, finally, of interest to understand the reason why the resonances here considered are 
broad and narrow relative to each other. Out of the five channels spanning the low-energy 
sector, the two 1 S+ singlets with spin configurations |0, > c |0, > n and |0, > |2,0 > n 
interact mostly with the 3 E+ triplet with spin configuration |1,T > |1,1 > n . In Fig. 3 the 
magnetic field dependence of the energy of the last singlet bound state (with 38 nodes) is 
plotted (dashed line) relative to the threshold of the triplet channel (dotted line) for the case 
when these channels are completely decoupled. These curves cross each other at (about) 550G. 
The triplet channel, in addition, is found to have a virtual state at threshold [24], as the large 
negative value of the background scattering length ab g indicates (cibg — —3790 a.u. when 
calculated with this channel only). [This value results from our ab initio calculation, where 
the only slightly adjusted parameters are the coefficients of the R~ 6 tail of the electronic 
energy curves. It differs by about 40% from that reported in the literature [23], where a 
larger number of adjustable parameters is however used. This difference in ab g does not affect 
appreciably the values in Table I but for the function ap(B).] 
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Out of the two singlet channels, it is always possible to find a linear combination that 
decouples from the triplet virtual state, thus crossing threshold at (about) 550G. The other 
combination which couples with the triplet is forced, correspondingly, to have an avoided 
crossing when the interaction between singlets and triplet is restored (full line). The combi- 
nation which decouples from the triplet results thus in a "narrow" resonance (with a small 
broadening provided by the coupling to the rest of the continuum states), and mantains its 
singlet character through the crossing at (about) 550G. The combination which couples with 
the triplet, on the other hand, acquires (almost) full triplet character past the avoided cross- 
ing. For large magnetic field, the energy of the latter tends asymptotically to the broadening 
(Ma^ g ) _1 of the virtual state (long-dashed line in the inset of Fig. 3), thus crossing threshold 
quite slowly. The slow convergence to the asymptote eventually accounts for the "broad" 
nature of the resonance vs the magnetic field. 

In conclusion, by extracting the relevant information from a coupled-channel calculation for 
the broad and narrow FF resonances of 6 Li, a single-channel model with a contact interaction 
proves adequate to describe the BCS-BEC crossover with the broad resonance. In this way, 
a connection is established between the BCS-BEC crossover with trapped Fermi atoms and 
analogous crossovers with other systems. Trapped Fermi atoms are thus ideally suited for 
studying the BCS-BEC crossover in a universal fashion, irrespective of mechanisms specific 
to these simple systems. 

* * * 
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